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ABSTRACT 

Observations of line ratios in the Milky Way's warm ionized medium (WIM) suggest that pho- 
toionization is not the only heating mechanism present. For the additional heating to explain the 
discrepancy it would have to have a weaker dependence on the gas density than the cooling rate, An^. 
Reynolds, Haffner & Tufte (1999) suggested turbulent dissipation or magnetic field reconnection as 
possible heatmg sources. We investigate here the viability of MHD-wave mediated cosmic ray heating 
as a supplemental heating source. This heating rate depends on the gas density only through its linear 

— 1/2 

dependence on the Alfven speed, which goes as ne . We show that, scaled to appropriate values of 
cosmic ray energy density, cosmic ray heating can be significant. Furthermore, this heating is stable 
to perturbations. These results should also apply to warm ionized gas in other galaxies. 



1. INTRODUCTION 

Observations of [S II] /Ha and [N II] /Ha line intensity ratios in the WIM show a spatial variation with distanc e 
from the galac tic midplane \z\ - larger line ratios are seen further from the disk (see Reynolds, Haffner fc Tufte] (|1999), 



Haffner et al. ( j2009) ). Such variation might be explained by variations in the ionization parameter U , the ratio of 
photon density to gas density. However, this would not explain the additional observation that the [S II]/[N II] ratio 
remains nearly constant with \z\. Under WIM conditions, variations in U inevitably produces larger changes in Sulphur 
(which can be either in the form of SII or SHI) compared to Nitrogen (which almost always appears as Nil), due to 
their different ionization potentials. 

These observations may be explained by a spatial variation of the electron temperature Tg. An increase in Tg with 
height above the disk could explain the enhanced [S II] /Ha and [N II] /Ha ratios. Also, because [S II] and [N II] have 
nearly the same excitation energy, the [S II]/[N II] ratio is nearly independent of Tg. So the near constant [S II] /[N II] 
ratio may also be explained this way. 

But how can this variation in Tg be explained? If only photoionization heating is important, then this increase in 
To can potentially be accommodated by hardening o f the spectr um away from the disk mid-plane. However, a hard 
spectrum is inconsistent with Hel A5876 observations ( Rand|1997[ |1998p . On the other hand, if there were a secondary 
heating mechanism with a weaker dependence on electron density than the dependence of photoionization heating, 
such heating would dominate far from the disk, where densities are low, and we would see a variation in that could 
explain the observed line ratio variations. 

Many such supplementary h eating mechanisms have been proposed, su ch as photoelectric heating f rom dust grains 
jWeingartner & Draine (2001)), magnetic reconnection (Raymond' (1992)), and turbulent dissipation (Minter&Span- 
per, ( ,1997p ). We study nere the possibility of cosmic ray heating. The process was outlined by Wentzel^ ( |197l[) , but 
not applied to the WIM (which had not been discovered at that time) . If a cosmic ray population has a bulk velocity 
faster than the l ocal Alfven spe e d VAi = B / \/AiTPi , ma gnetohydrodynamic Alfven waves are generated and exhibit 
unstable growth (Wentzel ( 1968 1 , [Kulsrud fc Pearce ( 1969) ). In a steady state, these waves are damped by some other 
process(es), transierrmg energy to the gas. In this way a cosmic ray density gradient can indirectly heat the plasma. 
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We will see that the resulting gas heating rate is proportional to rig , and is of the order required to exp lain the 
necessary tempera ture variations. Note that this process is quite different from collisional cosmic ray heating (Spitzer 
fc Tomasko (1968)). A strength of this mechanism is that the heating rate depends on only a few parameters which 
are either observable or can be estimated. We first consider the nature of CR trapping in the WIM in ^ before 
considering the CR heating rate and its local stability properties in S|3] 

2. ALFVEN WAVE EQUILIBRIUM 

For cosmic ray heating to be in place we must ensure that Alfven waves are present in the WIM with enough energy 
to scatter the cosmic rays. To do this we determine the wave damping, which in the WIM environment is due to 
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ion-neutral friction and non-linear Landau damping. We then balance this damping with the cosmic-ray-induced wave 
growth to obtain an equilibrium condition. This condition determines the power spectrum in the waves. We can then 
derive a mean free path for the cosmic rays and determine if they are well-trapped. We could also use the equilibrium 
damping rate to determine the heating of the gas, but as we will see the heating rate depends only on the characteristics 
of the cosmic ray population, provided they are well-trapped. 

2.1. Ion-Neutral Damping in Nearly Ionized Gas 



We follow Appendix C of Kulsrud & Pearce ( 1969 1 . These authors assumed the gas is nearly neutral; we assume it is 
almost fully ionized. We begin with the force equations for the neutral and charged components of the gas respectively. 
Assuming the transverse velocities of each component are in the form of an oscillator Vi = ^^e*''^^"'^*), we have 



PnUj'^Vn = -il^„iWp„(i;„ - Vi) 
PiJ^Vi = PiUsl - iVinU!pi{Vi - Vn) 



(1) 

(2) 



These equations are essentially those of two coupled, damped oscillators, one of which is driven with frequency 
ujk = kvAii the natural Alfven wave frequency. The neutral-ion collisions are treated as a drag force parameterized 
by the collision frequency Vni- The movement of the ions by the Alfven waves is impeded by the neutral particle 
population. This nudges the neutral component to follow behind the oscillating ions, removing energy from the Alfven 
waves to do so. 

We can rearrange these equations into convenient matrix form by using Vni/i^in — Pil Pn- 



{up- - Wfc + iixlVni)Vi ~ iuJVniVn = 
- iixlV^nVi + {up + iujlVin)Vn = 

Setting the determinant of this matrix to zero gives us the dispersion relation for w. 

up + iVni{^ + —)uP — Ul'^U! — iVniUjf. = 
Pi 



(3) 
(4) 

(5) 



We can solve this perturbatively in the two limits uj^ ^ Vni, 2> Vni- Let us further assume we are in the 
e = Pn/ Pi ^ 1 limit, to match the WEVQ Let's first consider the ujk <C Vni case. To leading order in uj^/vni we have 



iz/„i(l -I- e)wo - ivni'^l = 



UJQ 



= UJkVl + (- 



The first order equation is 



Wq -I- 2w„i(l + e)uJoU}i — ujIujq = 
ujley/l + e + 2w„i(l + e)uJkVl + ewi = 



UJl 



2Vnr 1 



So in this limit the waves are damped at a rate 



Li. 



2vni 1 + e 

In the short wave limit, ^ i/„i, which is the relevant limit for CR-generated waves, we have 

upQ - uP^ujQ = ^ Wo = Wfc 



2ujkUJi + ivnieuj^ = 



=> OJi = e = 

^ 2 2 



and so we have damping rate 



(6) 

(7) 

(8) 

(9) 

(10) 
(11) 



Note that this is the same damping rate derived in Kulsrud & Pearce (1969), even though we are in the opposite limit 
of a mostly ionized gas rather than a mostly neutral one. 
The ion-neutral collision frequency is 



n„(cru) = \nn{<Jv) = \eni{av) 

m„ 2 2 



(12) 



Here n„ is the density of the neutral component and {av) is the average rate of exchang e of velocity per particle fo r 
ion-neutral collisions. For temperatures around 10"' K, we have {av) « 10^* cm'^ s~^ (see Kulsrud & Cesarsky (1971), 



^ Note that our t is the inverse of the one used in 



Kulsrud &; Pearce 



119691 
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De Pontieu, Martens &: Hudson (2001)). We assume here that hydrogen is the dominant neutral species; up to 10% 
ot the hydrogen in t he WIM is thought to be neutral ( Haffner et al. ( [2009 )); He should be mostly neutral, but has a 
lower collision rate (De Pontieu, Martens & Hudson (20011). So in the short-wave, almost completely ionized limit, 
the ion-neutral damping rate is a tunction ot only the density of the neutral component 



--n„{<7v) 



(13) 



2.2. Non-Linear Landau Damping 

In some circumstances, non-linear Landau damping may be comparable to or dominate over ion-neutral damping. 
Non-linear Landau damping occurs when ions ride along the envelope of a beat wave formed by two interfering Alfven 
waves. Ions whose random motions are slightly slower than the speed of this envelope will take energy from the waves, 
damping them. Ions that are slightly faster will give energy to the waves, but for a thermal distribution we expect 
there to be more of the slower pa rticles, and so t he net effect is a wave damping. The strength of this damping depends 
on the strength of the waves as (Kulsrud (1978)) 



NLLD — 



(14) 



Here, Vi = y^ksT/m is the thermal speed of the ions. 



2.3. Cosmic Ray Instability 

A cosmic ray traveling along a magnetic field line with speed v and pitch angle cosine fi will interact with an Alfven 
wave with parallel wave number k^ under the resonance condition 



k. = 



»0 



(15) 



Here, is the cosmic ray's relativistic gyroradius, Qq = eBo/mc is its nonrelativistic gyrofrequency, and 7 is its 
Lorentz factor. In other words, a cosmic ray and an Alfven wave are resonant if the wave's wavelength is roughly equal 
to the dist ance t h e cosmic ray travels along the B-field in one gyration. 



Wentzel ( 1968 



Kulsrud & Pearce ( 1969 ) showed that a population of cosmic rays whose bulk velocity is faster than 
the Altven speed will spur unstable growth in th e waves. If we have such a distribution of CRs /(x, p, t), the resulting 
growth rate can be written, in the wave frame, (Skilling (1971)): 



Pgrowth (^2: 



2k, 



dl 



lip 



^ip 



(16) 



The above holds for Alfven waves propagating nearly parallel to the background magnetic field Bz. From here on we 
drop the z subscripts. The delta functions encode the resonance condition for CRs travelling in both directions. 

We can rewrite this expression in terms of the cosmic ray gradient along the background magnetic field In the 
absence of any sources or sinks, the cosmic ray transport equation is 



dt 



flV 



d_l 

dz 



d_ 



(l-M^) 



dl 



(17) 



The scattering frequency v is related to the energy density £ of resonant Alfven waves (Kulsrud & Pearce (1969): 

H,)-'pkm-l^{§)\ (18) 



This collision frequency is expected to be very large compared to the cosmic ray dynamical timescale (a condition we 
must che ck later for consistency), so we can expand / in inverse powers of i', / = /q + /i + /2 + .... To lowest order, 
eqn. (17 1 becomes 



= 



To first order we have 



If we integrate both sides over fi, 



d_ 
dfi 

djo 
dz 



fiv 



d_ [(1-A^^) 
dfi 



V{fi) 



dh 
dfi 

dh 
dfL _ 



= 



dh 

dfi 



V dfo 

V dz 



(19) 



(20) 



(21) 
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We can now eliminate ^ from eqn 
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rgrowth(^) 



Plugging in equation (18) for we have 



^~ k{5B)l 
Integrating the delta function over jj, gives 

rgrowth(^) 



a/ 

dz 



V dz 



+ S[fip + 



(22) 



S [fip 



k 



f-OO /" 1 

/ 2'Kp'^dp / (i/i(l — lJ'^)pv 

Jo J -I 



dl 
dz 



5{ fip- 



+ d{^p + 



k 



Sir'^mfloVA 



dpv — [p -p^) 
oz 



(23) 



where we have denoted — mflQ/k. 
To make eqn. (23 1 look a bit simpler, let us rewrite the integral in terms of a unitless factor of order unity A(k): 



Aik) = 

and let's define a CR length scale by 
The growth rate is therj^ 



Pk 



ncR 



dpvf{p){p^ - pi) = —cncRA{k) 



dpPf{p)A7T{p^~pl), 0<A{k)<l 



Pk 



dncR _ ncR 



CR 



r. 



growth 



(fc) = 



dz 

k{SBy ~ 



CR 



m 



(24) 
(25) 
(26) 



Without specifying a cosmic ray distribution / we cannot say anything about A{k). As an example, consider a 
power law in momentum, f(jx.,p,t) — C{x,t)p~°' , with some lower momentum cutoff pc and normalization: 



ncR(x,i) 



A'Kp^f{y.,p,t)dp 



in 
a — 3 



cp: 



3-Q 



/(x,p,i) = ncR(x,t) 



a — 3 f p 



Anpl \p, 



e{p-Pc) 



(27) 



Then by definition ( 24 1 we get 



A{k) 
(a -3) 



a — 3 



c -Jpk 



Pc 



If we take the relativistic limit /? w 1 and denote max(a;fc, I) ~ yk this becomes 



P Pk 

X = — Xk= — 

Pc Pc 



A{k) = (a-3) 



or 



3 — a 1 — a 

, Q-3 



Vk 



2 "-3 2' 
a — 1 



Aik) 



a-1 \ka 

a— 3 ( kc 



(28) 



1-1 I k J 

where kc = ttiVLq/pc is determined from the lower momentum cutoff of the spectrum. 

^ In principle the quantity A{k) could vary in space, but we ignore this possibility here. Alternatively we could adjust our definition of 
LcR to include this effect. 
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2.4. Equilibrium Power Spectrum 
Now that we have the total damping and growth rates of the Alfven waves we can enforce an equihbrium condition 



- growth 



NLLD 







(29) 



Inserting our expressions ( 14 ) , and ( 26 ) into eqn. ( 29 ) , 



2TrmfloVAC ncR 



A[k) — Tin ^ 



v,k 



5_B_ 
l3 



or, rearranging terms, 



X = 



l3 



We solve analytically for X 

'SB\' 



2 Tin 

TT kVi 



[-l + ^/l+7^], 



n 



c rj ncR 
2 VA LcR rii Tl„ 



A{k), 



(30) 



where r.i = Vi/flo is the thermal ion gyroradius. 

It is informative to determine the relative importance of each damping mechanism. We can do this by looking at 
the quantity TZ. If TZ is small, the linear term (ion-neutral damping) dominates and 

^ 1 c ncR iloA{k) 



SB 
l3 



^. 2vA Ui kLcRT^ri 
while if TZ is large, non-linear Landau damping dominates, and 

2 



7^< 1. 




e ncR A{k) 
k^rjL. 



1/2 



■7^> 1. 



TT VA rii 



CR 



The transition between these two limiting cases occurs at 7?. = 1, or 

-1 

A{k) « 2.8 X 10-4 ( , ~ 

^ ' V cm/s 



V.Olcm-'V V 



"-CR 



10- 



L. 



CR 



kpc 



e 

X)5 



(31) 



(32) 



(33) 



where we have used eqn. (13) and taken (av) ~ 10^® cm'^ s^^. We introduce here a set of convenient fiducial values 
that we will use throughout this paper. For the power law spectrum with a = 4.7 this gives us a transition wave 
number fc* of 

.012^ = 3.6x10-15 cm-i ( ^] ( 



.012fcc 



PcC 



fiGj VGcV/ 



(34) 



330 

lu da 



at the fiducial values in (33 For k ^ k* , ion-neutral damping domina tes and the wave power is given by (31 ). For 



k ^ k* , non-linear Landau damping is dominant and the wave power is ( 32 ) . 

We are now in a position to check whether the cosmic rays are self-trapped, i.e. whether their mean free path 
to scattering by self generated turbulen ce is small compared to their scale height. This is a necessary condition for 
applyin g th e heating theory derived in § |2.5| and[3l The mean free path A is related to the scattering frequency given 
in eqn. ( 18 ) by 

V 



(35) 



Using eqns. ( 18 1 and ( 30 ) in eqn. ( 35 1 yields a relatively compact expression for A 

'-i + \/r 



A = 



TZ 



(36) 



The mean free path in pc given by eqn. ( 36 ) is plotted as a function of p in units of the cutoff momentum pc in 
Figure [ij The figure spans the transition from Landau damping dominated at low momentum to ion-neutral friction 
dominated at high momentum. It appears from Figure [l]that cosmic rays of energy even several hundred times the 
cutoff energy arc quite well trapped (A ~ several pc). Most of the cosmic ray energy lies in the trans-relativistic ~GeV 
regime. In the ISM, the spectra turns over at p^ ~ 10 MeV due to Coulomb coo ling . Thus, the regime of interest is 
pIPc ~ 100. At the upper end of the trapped range, where eqn. (31 1 holds, eqn. (36) can be written in the form 



A 



Lm VA 



rii 



(37) 



-T 7S ^-^4.0x10-^ 

LcR TT ilo c ncRA{k) 

^ We note here that the value of fc* depends heavily on these quantities, particularly the gas density rii. In fact, for some values there 
is no region of fc-space where non-linear Landau damping dominates. 
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Figure 1. Cosmic ray mean free path to scattering by self generated turbulence, calculated from eqn. | |36[ l with the parameters set equal 
to the fiducial values. The mean free path is given in pc^aud the cosmic ray momentum is given in terms of the cutoff momentum pc- 

where in the last expression we have used eqn. ( |28[ ) and set all parameters to their fiducial values. 

Figu re [T ] can also be used to check that the waves are small amplitude and well described by linear theory. From 
eqns. (18) and (35l was can see that {6B/B)1 ~ r^/A. Cosmic rays of energy a few hundred GeV and less have 



gyroradn of order iOs of AU or less, showing that SB/B <C 1 even at low momenta where the mean free path is short. 

2.5. Heating Rate 

We can now determine the heating rate of the WIM due to the dissipation of Alfven waves created by cosmic ray 
streaming. To do this we want to integrate the time-derivative oi £{k) — SB'^/Snk over all wave numbers k. We know 
that the time-dependence of SB in an Alfven wave is 



SB oc e-"^* 

UJ = UJR + iPdamp 



and so 



£ik) 



SB?. 
8nk 



d£{k) 
dt 



(38) 



(39) 



(note that Fdamp < 0). 

Let us remove the assumption of a power law spectrum and a relativistic limit and go back to any general distribution 
/(x,p, t). Let us also remove any assu mpti ons about damping mechanisms, and only assume we have equilibrium for 



the Alfven waves. Then from equation (26) and Fgi-owtii = Tdamp, 



dfc2rda,„p(fc)-^ = / dfcFd 
Jo Jo 

mfloVAC d 



H 



dk- 



2fc2 



dz 



2Fdamp(fc)fc'^ LcR. 



(40) 



(41) 



In the last step we have rewritten A{k) and Lcr in terms of their original definitions (24) and (25). We reformulate 
this double integral with a change of variable from k to Pk, dpk — —mflQ/k^dk. Let us also write va = va^ such that 



dpk / dpAirfip) V (p) [p^ - pI ) 



(42) 



Finally, let's exchange the order of the integrals by recognizing that the double integral is over all {p,Pk) under the 
constraint < pfc < p. 

dpi dpkiTTf{p)v{p){p'^ -pI) 
'0 Jo 



H 



1 



-VA • V 



dpATrp^f{p)v{p) 



1 



dp4Trp^f{p)v{p) 



These integrals are now very simple - they correspond to the total cosmic ray pressure 

1 



^CR ^ 3 y dpAnp f{p)v{p) 
We therefore obtain the very simple expression for the cosmic ray heating: 



(43) 



(44) 



H = -VA • VPcR 



(45) 
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in agreement with Wentzel 



1971). 



The heating rate is simply the cosmic ray pressure gradient times the Alfven speed. Even with out the above 
calculation we know this must be the solution, since we require an equilibrium fo r the waves and (|45| is al ways the 
rate at which cosmic rays give energy to the Alfven waves regardless of damping (McKenzie & VoeE( 1981 )). So, as 
hinted at in section [2j we require only t hat the Alfven waves are in equilibrium arid the cosmic rays are well-trapped 
to know that the cosmic ray heating is (45). 

3. APPLICATION TO OUR GALAXY 
3.1. Observations 

Let us carry through the dependence of va on the magnetic field and ion density, and pick som e representativ e 
values. To write the cosmic ray pressure in terms of the energy density, we use Pqr — 0.45i?cR from Ferriere (2001 ). 
Then we can estimate the heating rate in the WIM: 



w 5.2 X 10 



-28 erg Ei 



CR 



B 



cm3 s eVcm fJ,G 



L 
kpc 



10- 



-1/2 



(46) 



If we assume the cosmic ray energy density and magnetic energy density fall off with the same scale height L, we can 
determine the dependence of this heating rate on height z from the galactic plane 



Biz) = Boe 



-\z\/2L 



H{z) « 5.2 X 10- 



-28 erg EcR,Q 

cm3 seVcm^^ mG V^pc 



Bo 



L 



10- 



-\z\/L 



1-3 



-1/2 



-3\z\/2L 



(47) 

(48) 



Let's compare th is to the heating that would be nec essary to explain the inferred temperature profile T{z). Following 
the prescription inlReynolds, Haffner & Tufte (1999) we find T{z) in our model by solving a heating-cooling balance 
equation 



Gonl + G^n-^/^ = Ar 



(49) 



where each term represents, from left to right, photoionization heating, cosmic ray heating, and the cooling rate. The 
temperature dependence of the electron density rig, the cooling function A, and the photoionization heating Gq are 



,(|z|) = 0.125T°-4V""-^e-l"l/i'^P= cm-3 
A = 3.0 X IQ-^^Tl^ erg cm^ s^^ 

Go = 1.2 X 10"24j,-0.8 gj.g ^^3 g-l 



(50) 
(51) 
(52) 

T4 denotes the temperature in units of 10** K, and / is a filling fraction describing the amount of ionized Hydrogen, 
which we will set to f{z) = Min[0.1el^l/^^°P'^, 1]. Note that our assumptions imply that VPc < pg at all z, and is 
consistent with hyd rost atic balance. 

We then solve eqn 49 to obtain the model profile T(z) and compare it to the profile infered from line ratio observations. 
We adjust G3 to fit tne model curve to the data, and we find, for L = 2 kpc: 

G3,fit « 1.2 X 10-27 ^^e-3l-l/4000 PC (53) 



=^ Hfit — G 



-1/2 _ 



= 1.2 X 10" 



cm^/2 
erg 



10- 



-1/2 



-3|2|/4000 pc 



(54) 



See figure [2] (blue points) for this fit. 

Comparing this to equation (48) with Lqr = 2 kpc, we see that cosmic ray heating is sufficient to explain the 
observed line ratios if the magnetic field and CR energy density normalizations are high enough: 



eVcm"^ pG 



46 



(55) 



Is this the case? The m agnetic field an d CR energy density in the solar neighborhoo d are about = 5 pG 
and £'cR,o = 1-8 ev cm~'^ (Ferriere (2001)). Beuermann, Kanbach & Berkhuijsen (1985) showed that synchrotron 
emissivity in the galactic spiral amis is auout 4 times greate r than m the mterarni regions. Synchrotron emissivity 
depends on the field and CR density as ( Pfrommer fc Enfilin] p004)) 



(56) 



for a power law CR density. We utilize our equipartition assumption in the last step. A spiral- arm enhancement of ji, 
by a factor of 4 therefore implies the i?-field increases by a factor of about 1.4 and the CR density increases by about 
1.9 (for a = 4.7). We might then expect the product of the i?-field and the CR density in the Perseus arm to be about 



Er 



CRfiBo ~ (1.8 eV cm"-')(5 /xG) * 1.4 * 1.9 w 24 iiG eV cm 



-3 



(57) 
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Figure 2. Temperature versus height for the WIM. The black line is derived from observations of line ratios in the Perseus spiral arm. 
The blue points show the fit solution using the above parameters. The red points show the profile using parameters based on observations. 



This faUs short of the requirement from equation ( |55[ ) but stiU comes close to reproducing the inferred temperature 
profile, as shown in the red points in figure [2] We also have not incorporated the orientation of the magnetic field. We 
have assumed a vertical field, but the actual field in the WIM may be much more random. This might be accounted 
for with an effective efficiency parameter. 

One complication these estimates do not fully take into account is the multiphase nature of the ISM: the WIM has 
a low 20%) filling factor, and above the disk is interspersed with hot diffuse coronal gas. This can cause local 

1/2 

variations in Alfven speed, and thus cosmic ray pressure. As long as B increases more slowly than Ue (no te that in 
the cooler d iffuse ISM sampled by HI lines, B as measured by Zeeman splitting does not scale with density (Crutcher 
et al.|[20T0 )), va will be reduced in the WIM relative to coronal gas, and thus Pc is higher. Thus, the relevant length 
scale tor cosmic ray pressure gradients could be the cloud size, rather than the global scale height we have adopted; 
this leads to larger heating ratef[^ On the other hand, if clouds are sufficiently small (Lcioud ^ A ^ lOpc), then the 
CRs smooth over these inhomogeneities and the global gradients are appropriate. Ha observations and photoionization 
modeling indicate that the WIM is likely to have both a s mooth and a clumpy component which fluctuates on a wide 



range of length scales, but there is no consensus picture (Haffner et al. 2009). We regard these issues as beyond the 



scope of this paper, but such considerations are illustrative of possible variations in the cosmic rate heating rate 



3.2. Stability 

We must check that the heating from the cosmic ray pressure gradient is stable under perturbations. If the heating 
increases compared to the cooling for a small element of gas perturbed to a higher temperature, the heating is unstable 
and we would get thermal runaway. To determine if this is the case, we need the change in the heating rate for a given 
perturbation. 

5(An2 - i?) = 5{Lon^T^-^ - Hon^T-^ "^ - HsBEcru-'^/'^) (58) 
In terms of perturbed quantities this becomes 
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(59) 



Making use of the fact that the initial state was in thermal equilibrium, Lt^ii?T^ = Han^T "-^ + H^BEcrti we 
can write this most generally as 



6{An^ - if) = ^ [l.9Lon2Ti-9 + O.SHon^T-"-^] 
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(60) 



* It also implies most heating occurs when the cosmic rays exit the cloud, when v\ ■ VPc < 0. 



Without specifying the perturbation this is as far as we can go. Once we relate the perturbed quantities with ST, we 
can determine whether the heating is stable or unstable. 

Let us consider an isobaric perturbation perpendicular to the magnetic field lines. The total pressure remains 
constant 

(61) 



^tot = Pg + Pb + Pgr = const. 



The above terms represent the gas pressure, magnetic pressure, and CR pressure respectively. The field lines are 
compressed along with the gas, so 

SB _ Sn 

"b" ^ V 



Let us further assume that the cosmic rays respond adiabatically: 

SEcB. _ SPcR _ Sn 

- — - 7CR — 

EcR Pgr. n 

Then we can use the constant pressure condition to relate Sn and ST. From Pg cx nT and Pb oc we get 

SPtot - SPg + SPb + SPcR = 



(62) 



(63) 
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(64) 



Putting this all together into equation ( 60 1 gives 
ST 



S{An' 



H) 



-1/2 



— [l.9Lori'Ti-9 + 0.8ffon'r-°'- H^BEcRn 7CR 



Pg + 2Pb + 7crP( 



CR 



(65) 



If the term in the brackets is positive, a small increase in temperature causes the change in cooling to outweigh the 
change in heating and the perturbation is stable. If the term in brackets is negative, it is unstable. But note that 
both parenthesised factors in the third term must each be less than one. Also, by the thermal equilibrium condition, 
H^BEc^nT^/^ must be less than L^n^T^-^ . The first term must therefore be of higher magnitude than the third term, 
and so the expression in the brackets is positive and the heating is stable. 

This is perhaps easier seen if we denote the total cooling by C, and the fraction of the total heating due to 
photoelectric heating by < x < 1. Then, 



S{W -H) = C 



ST 
~T 



1.9 + 0.8X- {l-x) 



IGR 



Pn 



Pg + 2Pb + icrPgr 



(66) 



Now let us consider an acoustic perturbation. Equations (60), (62), and (63) still hold. Pressure is no longer fixed, 
but the gas responds almost adiabatically. As such 
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(67) 
(68) 



We therefore have 



5(An2 -H) = 



ST 
~T 



1.9Lo?^ T 



Q.m^n^T 



HsBEcRn 



-1/2 3/2 - 7CR 



79-1 



or, in terms of C and x, 



S{W -H) = C 



ST 
~T 
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7,-1 



(69) 



(70) 



By a similar argument as before, this is also stable. 
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4. SUMMARY AND CONCLUSIONS 

The gaseous disk of the Milky Way has a warm ionized component (WIM) with scale height several times that of 
the predominantly neutral component. The magnetic field and cosmic ray components have similar thickness. Thick 
layers of warm ionized gas, and extended nonthermal emission, are seen in other galaxies as well ( :Haffner et al.^2009( ). 

It is widely accepted that starlight photoionizes and heats the WIM. Nevertheless, there is evidence tor a supplemental 
heating mechanism. Detailed reconstruction of the WIM vertical temperature profile in the region of the Perseu s spiral 
arm shows an increase in temperature with height that cannot be explained by radiative heating alone (Reynolds, 
Haffner & Tufte (1999)). These authors showed that these observations can be explained by an additional heatmg 
mechanism w ith a weaker d ensity dependence than the n ? dependence of radiative heating. Heating by magnetic 



recon nection (|Raymond| (|l99 2)), dissipation of turbulence (Minter & Spangler (1997)), and photoelectric heating by 
dust ( Weingartner &: Uraine| (j2001j)) have all been invoked. All three are feasible on energetic grounds, but the rates 
of the first two, m particular, depend on many unknown factors and are quite uncertain. 

In this paper, we estimated the heating rate due to dissipation of waves excited by streaming cosmic rays. When 
the cosmic rays are well scattered by this self-generated turbulence, the heating rate depends only on the cosmi c ra y 
pressure gradient projected along the local magnetic field direction and th e magnitude of t he Alfven speed (eqn 



Cosmic ray heating of the interst ellar medium was discussed in ge neral in Wentzel ( 1971 ) and is incl uded in mooe 
of supernova driven shock waves (fVoelk, Drury fc McKenzie (1984)), cosmic ray driven galactic wind s ( Br eitschwerdt 



45) 



McKenzie k Voelk| ( | 1991 [ ), [Everett et al.| ( |2dD5 )), and dittuse mterstellar clouds ( [Everett fc Zweibelj ( |20lip ) but up to 
now does not appear to have been considered tor the WIM. In §2 we showed that when wave excitation by streaming 
is balanced against nonlinear Landau damping and ion-neutral friction, the resulting wave amplitude, while small 
enough to allow treating the waves in the small amplitude approximation, is large enough to scatter the majority of 
cos mic rays many times over one pressure scale height (Figure 1). This justifies the frequent scattering limit we used 
in j |3.1[ to estimate the cosmic ray heating rate for the WIM (eqns 46 and 48 ) and show that adding it to the thermal 
equilibrium model of the WIM for Perseus Arm conditions produces a reasonably good fit to the observations (Figure 



57|), the height dependence - which 
leads to a temperature 



2). Although the heating rate coefficient is about a factor of 2 too small (eqn. 
follows from the height dependence of the magnetic field, gas density, and cosmic ray pressure ■ 
vs height relation of the correct shape. We regard this, and matching the inferred size of the supplemental heating 
rate to within a factor of two - as confirmation that cosmic ray heating is a viab le supplementary heat source for the 
WIM. Cosmic ray heating also seems to be a thermally stable mechanism (3.2), at least under the assumptions we 
considered. 

The results in this paper should be generally applicable to warm ionized gas in other galaxies. In cases where 
synchrotron emission is detected or other estimates of the cosmic ray and magnetic field energy densities are available, 
it should be possible to estimate the magnitude of cosmic ray heating. It is important that the gas be diffuse and that 
the ionization fraction be high; in weakl y ionized clouds, for exam ple, ion-neutral friction is so strong that the cosmic 
rays are not well coupled to the medium (Everett & Zweibel ( 2011[ )). And, as long as the cosmic rays are well scattered, 
their pressure gradient along the ambient magnetic field exerts a force which may be important in determining the 
scale height of the gas and in driving an outflow even when the thermal speed of the gas is well below what is needed 
for escape. 
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